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ABSTRACT 

We studied the distribution of dust clouds around the central engine of NGC 
1068 based on shifted-and-added 8.8 - 12.3 /im (MIR) multi-filter images and 3.0 
- 3.9 /im (Z-band) spectra obtained with the Subaru Telescope. In a region of 
100 pc (1'.'4) around the central peak, we successfully constructed maps of color 
temperatures and emissivities of the MIR and L-band continua as well as the 
9.7 /xm and 3.4 /xm dust features with spatial resolutions of 26 pc (0'.'37) in the 
MIR and 22 pc (0'.'3) in the i^band. Our main results are: 1) color temperature 
of the MIR continuum scatters around the thermal equilibrium temperature with 
the central engine as the heat source while that of the L-band continuum is higher 
and independent upon distance from the central engine; 2) the peak of the 9.7 /im 
silicate absorption feature is shifted to a longer wavelength at some locations; 3) 
the ratio of the optical depths of the dust features is different from the Galactic 
values and show complicated spatial distribution; and 4) there is a pie shaped 
warm dust cloud as an enhancement in the emissivity of the MIR continuum 
extending about 50 pc to the north from the central engine. We speculate that 
material falls into the central engine through this cloud. 

Subject headings: galaxies: active — galaxies: individual (NGC 1068) — galax- 
ies: nuclei — galaxies: Seyfert — infrared: galaxies 
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1. INTRODUCTION 

Physics around black holes, including accretion of material and launching of jets, has 
been a fundamental question in astronomy for a long time (e.g. Eardley & Press 1975). It is 
not yet clearly known how accreting material disposes of its angular momentum while falling 
into a black hole, and how and where a jet is launched from around a black hole. With a jet 
of 100 pc scale observed (e.g. Gallimore et al. 1996), active galactic nuclei (AGNs) of Seyfert 
galaxies are the best targets to attack the mystery of accreting material and jets. 

There are two classes of Seyfert galaxies; Seyfert 1 and Seyfert 2. They have similarities 
and differences in, for example, their visible spectra, which are commonly known to be 
explained with the unified model (Antonucci 1993). According to the unified model, an 
AGN has a black hole at its heart. Material falls down into the black hole through an 
accretion disk while emitting X-ray and UV light. Throughout this paper, we define the 
system of the black hole and the accretion disk as "central engine" following Galliano et al. 

(2003) . The unified model proposes three classes of clouds around the central engine. First, 
narrow line regions (NLRs) spread over ~ 100 pc and emit permitted and forbidden lines 
with velocity widths of < 1000 km sec~^. Second, broad line regions (BLRs) in the vicinity 
of the central engine emit permitted hnes of ~ 10000 km sec~^ velocity widths. Third, an 
opaque compact but geometrically thick dusty torus around the central engine hides the 
direct view to the central engine and their BLRs from our line of sight when the central 
engine is viewed through the torus. The dusty torus is the key to the unified model because 
it is the one which distinguishes Seyfert 1 and Seyfert 2: Seyfert 2 nuclei are viewed through 
the dusty torus and BLRs are only seen in scattered light. 

NGC 1068 is one of the nearest Seyfert 2 galaxies. Many observations have revealed 
various features that confirm the unified model; an accretion disk of 1 pc scale (Gallimore 
et al. 1997, 2001), the Seyfert 1 nucleus seen in polarized light (Antonucci & Miller 1985), 
and NLRs (Evans et al. 1991). However, no clear image of the dusty torus has been obtained 
(e.g. Bock et al. 2000; Tomono et al. 2001; Rouan et al. 2004). Only recently, Jaffe et al. 

(2004) found a 2.1 pc scale emission component of warm (320 K) dust surrounding a more 
compact hot (> 800 K) dust component with the Mid-Infrared Interferometric Instrument 
coupled to the Very Large Telescope Interferometer. They identified the hot component 
as the inner wall of the dusty torus. It is a step towards clear images of the dusty torus 
including its temperature structure, density distribution, and relations with structures in 
larger and smaller scales. 



^ Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa- 
tory of Japan. 
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On the other hand, the inner radii of the dusty tori are being measured with reverbera- 
tion mapping. This method measures lag time of hght curves of an AGN in the F-band and 
that in the i^-band. In reverberation mapping, the measured lag time is interpreted as the 
time duration in which light from the AGN reaches the inner wall of the dusty torus. Multi- 
plying the lag with the speed of light, the inner radius of the dusty torus Rin can be obtained. 
Minezaki et al. (2004) found a good correlation between lags and absolute magnitudes in the 
F-band in a number of Seyfert 1 galaxies. It implies that Rin is governed by luminosity of 
the central engine. Adopting their results, minimum for NGC 1068 is estimated to be 
~ 0.045 pc from my ~ 11.5 mag (Sandage 1973). The estimated radius is consistent with 
the upper limit of 1 pc measured by Jaffe et al. (2004). The method seems to be promising 
in unveiling the innermost structures of dusty clouds around AGNs. However, we have to 
note that the results do not necessarily require the dusty clouds around the AGNs to be in 
the shape of a torus. It is still critical to obtain direct images of dusty tori. 

It is best to look for a dusty torus in the infrared (IR). First, from thermal equilibrium, 
it is expected that dust in a torus of 10 pc scale radiates mainly around 10 ^m. (Pier & 
Krolik 1992). Next, there are absorption/emission features of carbonaceous dust at 3.4 /im 
and of silicate dust at 9.7 /xm. Comparison of the dust features may be used to probe the 
temperature structure of a dusty torus as has been suggested by Imanishi (2000). Because 
i^band emission is dominated by dust at ~ 1000 K, the 3.4 /xm absorption feature should 
originate from the dust clouds in a few pc from the central engine (Pier & Krohk 1992). 
On the other hand, the 9.7 /xm absorption feature only measures the extinction towards the 
outer region of the dusty torus because dust of ~ 300 K, which is the dominant source of the 
mid-infrared (MIR) emission, is located further out from the heating source. Therefore, the 
combination of MIR spectra (8 - 13 /xm) and Z-band spectra (2.8 - 4.2 ^m) would probe 
the temperature distribution of dust in the line of sight through the dusty torus. 

This paper reports the results from observations of NGG 1068 in the MIR and in the 
i^band with the Subaru Telescope, which is best suited for this study with the superb image 
quality in IR. Two instruments, the Mid-Infrared Test Observation System (MIRTOS) and 
the Infrared Camera and Spectrograph (IRCS), were used with the fine pixel scales optimized 
for the image quality of the telescope. The remainder of this paper is organized as follows: 
Sections 2 and 3 describe observations and data reduction in the MIR and in the i^band, 
respectively; Section 4 discusses emission mechanisms of the continua, spectral and spatial 
change of the dust features, and distributions of warm dust clouds; and § 5 summarizes the 
results. Throughout this paper, we assume the redshift of 0.0038 for NGC 1068 (Bottinelli 
et al. 1990). One arcsecond corresponds to 72 pc at the distance of NGC 1068, assuming 
the Hubble constant of Hq = 75 km sec~^ Mpc~^. Galliano et al. (2003) measured positions 
of the various spatial structures in the NGC 1068 nucleus observed at different wavelengths 
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and found that the central engine is at the IR peak. We employ their registration to compare 
the results in the MIR and in the L-band as well as those in the hterature. 



2. MIR DATA 
2.1. Observations 

MIR imaging observations of NGC 1068 with the "sihcate filters" (A/AA ~ 10) at 7.7, 
8.8, 9.7, 10.4, 11.7, and 12.3 fim were made on 1999 December 31, 2000 January 9, and 2000 
January 18 UT. We used the MIRTOS (Tomono 2000; Tomono et al. 2000) mounted at the 
Cassegrain focus of the 8.2-m Subaru Telescope. Pixel scale was 0'.'067, corresponding to a 
projected distance on the sky of 4.8 pc at the distance of NGC 1068. Sequences of 96 images 
of 31 msec integration were obtained. Between the sequences, the telescope was nodded 
typically 30' to the north to acquire background images. In this study, the data described in 
Tomono et al. (2001) were reduced again with the higher computing power available today. 

2.2. Shift-and-add and registration of images 

Many short-exposure images with 31 msec integration time were shifted-and-added to 
construct high-resolution images in each filter. First, background images, which were made 
with averaging the short-exposure images in the nodding pair, were subtracted from the raw 
images. Second, bad pixels and detector array specific features were corrected. Third, the 
peak position of the object in each image was measured with a precision of one tenth of a 
pixel with Gaussian fitting. The FWHM of the Gaussian was fixed to be the same as that 
of the Airy pattern to pick up the brightest speckle. Fourth, each pixel of the images was 
divided to 10 times 10 pixels resulting with a pixel scale of 0'.'0067. Fifth, each image of 
the smaller pixel scale was shifted by registering the peak position measured in the third 
step. Finally, the shifted images were averaged after some images with low peak values were 
removed. Table 2.2 summarizes total integration time of the short-exposure images used and 
the ratio of the images used to those acquired. 

The fiux was calibrated with the standard star a Ari (Cohen et al. 1999). Flux uncer- 
tainty was estimated at about 9% from the measurement of the other standard star a CMa 
(Cohen et al. 1995). 

The shifted-and-added image of each filter was registered and convolved to have the 
same FWHM of ~ 0'/37 (26 pc) as follows. The peak position of the object for each filter was 
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measured by Gaussian fitting of the pixels with flux more than 75% of the peak. The images 
were registered assuming that the peaks for each filters are at the same location on the sky 
(Galliano et al. 2003). The images are then Gaussian convolved so that all the images have a 
similar FWHM. Figure 1 shows the spatially integrated SED within a 3'.'8 diameter aperture 
on the registered narrow-band images. The ISO-SWS spectrum in the 14" x 27" aperture 
(Lutz et al. 2000; Sturm et al. 2000) is also shown for comparison. The flux density measured 
from the MIRTOS images are consistent with the flux density of the continuum in the SWS 
spectrum. 



Wc fitted the MIR spectral energy distribution (SED) with a grey body emission model 
including the silicate feature to estimate the temperature and amount of MIR emitting dust, 
as well as optical depth of the silicate absorption/emission feature. It should be noted that 
the extracted parameters do not represent the entire physical conditions but only characterize 
distributions of dust in some aspects as described below. The observed SED was fitted with 
a model flux 



where Bx{T,X) is the Planck function of temperature T and A;(A) is the optical depth of 
the silicate feature obtained from the IR excess towards Cep (Roche & Aitken 1984) and 
normalized at 9.7 /xm. 

The /i Cep extinction curve is used in this work to compare the silicate feature with the 
3.4 /xm feature as has been done by Roche & Aitken (1984). The 7.7 fim image was not used 
because the IR excess data does not cover the short wavelength. The peak wavelength of 
the silicate absorption feature seen in the SEDs in Figure 1 is different from that observed 
towards /i Cep: the peak wavelength of the absorption is shifted from 9.7 /im to 10.4 /im 
especially at the central peak. Because of this, the 10.4 iim image was not used in the SED 
fitting. See § 4.2 for a discussion on the change of the peak wavelength. 

The variables Tcont, ^cont, and T9.7 were treated as free parameters for the fitting. Color 
temperature Tcont represents the temperature of warm dust grains emitting the MIR con- 
tinuum (see § 4.1.1 for a discussion). The factor Scont is a product of dust emissivity at 
10 /xm and beam filling factor of warm dust grains of Tcont- Throughout this paper, we call 
Econt as "emissivity". The mass of the warm dust grains per unit area can be estimated as 
Mdust — ^cont/Kabs, where Kabs is the absorption cross section per unit mass of dust grains. 



2.3. Grey body fitting 




(1) 
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Fig. 1. — SED and results of grey body fitting in the MIR (lower data points, curves and 
left axis) and fiux densities in a 3'.'8 aperture (upper data points and right axis) compared 
with the spectrum measured with the ISO-SWS (Sturm et al. 2000, right axis). Each spatial 
element is measured in an area O'.'IS long (P.A. 3°) by O'.'S wide corresponding to the L- 
band slitlet. Each data point is shown with the filter path band as a horizontal error bar 
and fiux uncertainty of la, including uncertainty in the fiux conversion factor, as a vertical 
error bar. The thick curves display the resulting grey body models (eq. [1]) on respective 
spatial elements. The dotted curve shows the underlying model continuum for the peak. 
The 10.4 //m data are not used in fitting (see § 4.2). 
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The dust model (Weingartner & Draine 2001; Draine 2003) with Ry — 5.5 tabulates the 
absorption cross section at 10 fim as Kabs — 0.40 cm^ g~^. Although the dust cloud might 
be clumpy, dumpiness does not affect the mass estimation when Scont < 1 and a homoge- 
neous temperature is assumed of warm dust that is emitting continuum. The other free 
parameter, rg.r, is the optical depth of the silicate feature at 9.7 fim towards the warm dust 
that is emitting continuum. A negative optical depth in our model represents emission of 
the silicate feature. Although the spectral index of dust emissivity n is usually uncertain 
ranging from 1 to 2, we assumed n — 1.6 following Cameron et al. (1993). Figure 2 shows 
spatial distributions of the parameters fitted on each pixel. The images were successfully 
fitted over a region of 50 pc to the north and 40 pc to the south from the central peak: un- 
certainty in fitted parameters are Arg.y < 1, ATcont < Tcont/'^, and Aecont < £cont/2. Using 
the secondary reference star a CMa rather than a Ari, resulting Tcont is within 1.5 K, rg 7 
becomes about 0.1 deeper, and Scont becomes 10% stronger. To improve the signal-to-noise 
ratio of the parameters in the area surrounding the central peak, the data were also fitted 
after being convolved to have 1.4 times wider FWHM. 

Maiolino et al. (2001) pointed out that ground based measurements of the silicate ab- 
sorption feature at 9.7 fim might suffer from contamination in continuum with the uniden- 
tified infrared band (UIB) emission at 7.7 /xm. However, Sturm et al. (2000) detected little 
emission of 7.7 //m UIB towards NGC 1068 with ISaSWS in the 14" x 27" aperture. The 
3.3 fjLm UIB feature is not seen in the spectrum obtained towards NGC 1068 by Imanishi 
et al. (1997) with the 3'.'8 aperture nor in the L-band spectra in this work. In fact, the UIBs 
detected with ISOCAM (Le Floc'h et al. 2001) originate almost exclusively from the star- 
burst regions, peaking at a distance of 1 kpc (14") from the central engine. They indicated 
that the AGN contributes less than ~5% to the total integrated UIB emission. Therefore, 
we concluded that the contamination by the UIBs to the continuum is negligible in our field 
of view of of ~200 pc (2'.'7) around the central engine of NGC 1068. 

High temperature areas seen 0"34 from the central engine on east and west sides coincide 
with the second and third diffraction rings (0'.'29 and 0"41 radius) of the Airy pattern at 
8.8 nm. Between these diffraction rings, there are also diffraction rings at longer wavelengths 
which might affect the fitting. Nevertheless, a fitting on convolved Airy images does not 
reproduce a peak in fitted temperature. Therefore, we concluded that the high temperature 
region is not an artifact induced by the Airy rings. 

Flux was also spatially integrated over the area of each L-band slitlet (see § 3.2) and 
the integrated flux is fitted with the same model as described in eq. [1]. With comparable 
spatial resolutions achieved for the MIR shifted- and- added images and the i^band seeing 
limited spectra, no further convolution or deconvolution was performed. Results are shown 
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Fig. 2. — MIR image and grey body parameters (eq. [1]). The fit for the outer region is 
performed on images convolved with a Gaussian of FWHM same as the original PSF (the 
circle in panel a). Locations of the 5 GHz sources (Gallimore et al. 1996, open circles) and the 
[O III] clouds (Evans et al. 1991, filled circles) are indicated, (a) 11.7 //m image shown with 
5 contours in dex. The straight lines show positions and widths of the three i^band slits. 
(6) Tg.y map. Solid contours show rg.y between 0.1 and 1.5 (absorption; darker grey scale) 
with intervals of 0.1 while dashed contours show rg.y of and —0.1 (light), (c) Tcont(MIR) 
shown with contours between 160 K (dark) and 280 K (light) with intervals of 10 K. (d) Scont 
shown with five contours in each dex between lO"^ "^ (dark) and 10~^'^ (light). 
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in Figure 1. 

3. Z-BAND DATA 
3.1. Observation 

A seeing limited spectroscopic observation of NGC 1068 in the i^band was conducted 
on 2000 September 24 UT. We used the IRCS (Tokunaga et ah 1998; Kobayashi et al. 2000) 
mounted at the Cassegrain focus of the Subaru Telescope. The Z-band grism for 2.8 - 4.2 /im 
was used with the 0'.'3 slit, yielding the spectral resolution A/AA ~ 400. The length of the 
slit was 18". Pixel scale was 0'.'058 pixel"^. NGC 1068 was observed with the slit position 
angle (P-A.) of 3° from the north, with slit positioned at three positions to map the eastern 
area of the central peak: pointing offsets of ±6" along the slit and 0'/28 and 0'.'56 to the east 
were applied. Positions of the slits are illustrated in Figure 2 a. We call the sht at the central 
peak as "C", and the slits in the east as "0.28E" and "0.56E" according to the offsets to the 
east in arcseconds. Eight exposures of 2 seconds on-source integration time were obtained for 
each slit position. The sky condition was photometric throughout the observation. Airmass 
to the target was 1.09 ± 0.01. From the size of the image along the slit, spatial resolution 
was estimated to be O'.'Sl or 22 pc in FWHM. As a telluric standard, a spectrum of BS 813 
was obtained at the similar airmass (1.08). We selected a F-type star as a telluric standard 
because it is known that no strong stellar absorption hues are present in L-band for F-type 
stars. For flat fielding, spectra of dawn sky were obtained at the end of the night. 

3.2. Data Reduction 

The L-band data were reduced following the standard procedure: the 2-dimensional 
spectrograms were background subtracted by the dithering pairs, fiat fielded, and corrected 
for bad pixels using IRAF^. Subsequently, positional shifts of the spectrograms of each sht 
position were measured with cross correlations and confirmed to be not more than 0.2 pixels 
(O'/Ol or 3.2 A) in both the spatial and dispersion directions. The spectrograms were shifted 
in the dispersion direction on the 2-dimensional spectrograms to cancel wavelength shifts 
between the slit positions. Flux of the spectrograms were then calibrated with the spectrum 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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of the standard star BS 813, which was assumed to have an effective temperature of 7244 K 
(Johnson & Wright 1983). Wavelengths of the spectra were cahbrated using cross-correlation 
with a model atmospheric transmittance (Rothman et al. 1998). Then, the spectra in slitlets 
of O'.'IS along the slit were extracted to examine spatial variations. Finally, the spectra were 
convolved with a Gaussian with the FWHM of 8.2A in the dispersion direction to match the 
spectral resolution to that of the instrument and improve the signal-to-noise ratio. 

Figure 3 shows spectra of some slitlets in the slit C. It should be noted that the slopes 
of the continuum at A < 3.0 fim and A > 3.9 fim are sensitive to the atmospheric condi- 
tions during the observation. In the following analysis, we did not use the data in these 
wavelength ranges to avoid the uncertainty of the continuum slope. The periodic feature at 
3.46 - 3.56 /im in the spectrum of the central peak in Figure 3 is from the periodic telluric 
absorption in this wavelength range. Because the instrumental line profile changes upon 
spatial distribution of the illumination within the slit, the difference of the instrumental line 
profiles of the target and the standard star causes the residual features when the two spectra 
are divided each other (see the detailed discussion in Goto 2003). 



3.3. The 3.4 /^m absorption feature and the L-band continuum 

The absorption feature around 3.4 /im by carbonaceous dust is clearly seen in the 
spectra in Figure 3. Following Sandford et al. (1991), optical depth of the dust feature T3.4 
was measured as follows: 1) the continuum level at the absorption feature was estimated by 
connecting the continuum levels at rest wavelengths 3.23 and 3.64 fim with straight line and 
2) r3 4 was derived by averaging the optical depths at two sub peaks of the dust feature at rest 
wavelengths 3.38 fim (2955 cm~^) and 3.42 fim (2925 cm~^). The 3.48 fim (2870 cm^^) sub 
peak was not included in the measurements because of the residual telluric feature described 
in § 3.2. Figure 4 shows spatial change of t^a along with rg.y and rio.4. 

Color temperature Tcont and emissivity Scont in the i^band were estimated by fitting a 
grey body model 

Fx{X) = Scont ^ j^Q Bx{Tcont,^) (2) 

to the spectra at 3.00 - 3.23 //m and 3.64 - 3.90 //m. These wavelength ranges were chosen 
to avoid the 3.4 /im dust feature and the unstable continuum level at the both ends of 
the wavelength range due to changes of the atmospheric conditions. We assumed a dust 
emissivity index of n = 1.6 for comparison of the parameters with those in the MIR (§ 2.3). 
Spatial change of the color temperature in the i^band along with those in MIR are shown 
in Figures 5. 
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Fig. 3. — Z^band spectra on slitlets of (//IS long and C/S wide. Two vertical dotted lines 
at 3.23 iim (2955 cm~^) and 3.64 /xm (2925 cm~^) show the location of the sub-peaks of 
the 3.4 /im absorption features. The dashed lines connecting continuum levels at 3.23 and 
3.64 iim show continua for measuring T3.4 (Sandford et al. 1991). The dotted curves show 
model continuum (eq. [2]) with dust emissivity index oi n — 1.6, fitted for the data between 
3.0 and 3.9 pLm excluding the wavelength range of the dust feature (3.23 and 3.64 //m). The 
crosses in the left show typical la statistical uncertainties (vertical length) of each data point 
and the spectral resolution (horizontal length). Periodic features between 3.46 and 3.56 //m 
seen in some spectra are artifacts due to telluric absorption (Goto 2003). 
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Fig. 4. — Spatial variation of optical depths of the dust features for three slit positions. 
Open circles, filled circles, grey filled squares show T3.4, T9.7, and T10.4, respectively. Positive 
values of r mean absorption while negative values mean emission. Error bars are the la 
uncertainty. Error bars for T10.4, which are similar to that for rg.y, are not shown. The 
vertical line at bottom left shows the spatial resolution (FWHM) of the data. 
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Fig. 5. — Color temperatures of the continua in the i^band (open circles and upper abscissa) 
and in the MIR (filled circles and lower abscissa) on the three sht positions. Error bars show 
the la uncertainty of each measurement, which are less than the size of the symbols for the 
L-band measurements. The vertical line at bottom left shows the spatial resolution (FWHM) 
of the data. 
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RESULTS AND DISCUSSION 



4.1. 



Emission mechanisms of MIR and L-band continua 



Continua are detected both in the MIR and in the i^band over a region of 100 pc to the 
north and to the south from the central engine. Because Gratadour et al. (2003) detected 
no evidence of stellar activity in the vicinity of the central engine in their K-hand spectra, 
stellar activity should be negligible as a source of the extended IR continua. In the following, 
we will show that the MIR continuum is emitted from dust in thermal equilibrium with the 
UV radiation from the central engine while the continuum in the i^band is emitted by Very 
Small Grains (VSGs). 



Assuming a dust emissivity index n 1.6, Cameron et al. (1993) estimated dust temper- 
ature T in thermal equilibrium with a UV radiation field as 



where L is the luminosity of the UV source and R is the linear distance from the heating 
source. Figure 6 shows color temperature Tcont measured in the MIR and the estimated 
temperature for L = 2.2 x 10^^ Lq (Telesco & Harper 1980, at luminosity distance of 15 Mpc) 
as a function of the distance from the central engine. The distance is assumed to be the same 
as the projected distance on the sky. The vertical dotted fines in the Figure show spatial 
resolutions of data in the MIR and in the i^band, which are different because of different 
wavelengths and observation conditions. Because the light from the central engine, of which 
color temperature might be different from that for the extended continua, might contaminate 
the continua, the color temperatures measured around the central engine within the spatial 
resolutions are not accurate. 

In Figure 6, the measured Tcont in the MIR scatters around the estimated equifibrium 
temperature. This implies that, as a whole, the MIR emitting dust is mainly heated by the 
central engine and is in thermal equilibrium as has been suggested by Galliano et al. (2005). 



We interpret that the extended i^band continuum is emitted by VSGs among the fol- 
lowing three possibilities. First, the continuum is not dominated by the scattered light of 
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Fig. 6. — Color temperatures of continua measured on each slitlct in MIR (filled circles) and 
in the Z-band (empty circles). Horizontal axis shows projected distance on the sky from the 
central engine. Error bars show the la uncertainty of each measurement, almost of which are 
less than the size of the symbols for the Z-band measurements. The solid line shows thermal 
equilibrium temperature of dust (eq. [3]) heated by a UV source of L = 2.2 x lO-*^^ Lq 
(Telesco & Harper 1980). The spatial resolutions of the data and contamination of fight 
from the peak limit the color temperatures to be accurate only to the right of the vertical 
dotted lines (see text) . 
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the continuum from the central engine. Packham et al. (1997) measured spatially resolved 
polarization degrees in the J, and /^'-bands around the central engine and found that 
they are not more than than 5%. Second, the i^band continuum is not emitted from dust 
in thermal equilibrium. As shown in Figure 6, color temperature of the L-band continuum 
is higher than expected for dust in thermal equilibrium. Third, color temperature in the 
L-band remains roughly constant as far as 100 pc from the central engine; suggesting that 
hot grains emitting i^band continuum exist in this area even if their column density is not as 
high as those emitting MIR continuum. As Rouan et al. (2004) and Gratadour et al. (2005) 
suggested, the uniformity of the color temperature can be explained as emission from VSGs, 
which emit continuum with color temperature higher than their equilibrium temperature 
due to thermal fluctuation upon absorption of UV photons (Sellgren 1984; Sellgren et al. 
1985). 

If the continuum is emitted by VSGs, the incident UV photon flux does not determine 
color temperature but determines brightness, or emissivity in our grey body model, of the 
continuum. We actually found the change of emissivity of the L-band continuum upon 
incident UV photon flux from our grey body flttings, suggesting that the Z-band continuum 
is indeed emitted by VSGs. Figure 7 shows correlation between econt(-^band)/£cont(MlR) 
and Tcont(MIR). Thermal equilibrium temperature Tco„j(MIR) should correlate with the UV 
photon flux. On the other hand, since Scont is proportional to column density of dust grains 
emitting the continuum, £cont(-^band)/£cont(MIR) gives the ratio of dust grains excited to 
be VSGs to those in thermal equilibrium emitting continuum in the MIR. In more detail, 
energy F of UV photons absorbed by dust grains is balanced with the grey body emission 
with temperature T as follows: 

F^7^ j Scout (y^) ^a(T, A) d\ OC £contr^+", (4) 

where Econt is the emissivity of the dust cloud at 10 /xm. Division of two equations for the 
L-band and MIR continua leads to the following equation: 



Scont (-^band) F(L-band) 

£eont(MIR) " F(MIR) 



T,„„,(MIR) ]^+" 



7;ont(L-band) 



The solid line in Figure 7 shows the fltting result when assuming n — l.Q and adjusting the 
factor Tcont (L-band) x [F(MIR)/F(i^band)]^/^^+"') as a free parameter. Relation between 
the measured £cont (-^band) /econt (MIR) ^-nd the measured rcont(MIR) in Figure 7 roughly 
matches with the assumption in the range of almost two orders of econt(-^band)/£eo„t(MIR). 
The fltted factor (651 K) and the observed Tco„t( iv-band) ~ 560 K leads to F(L-band) / F(MIR) 
0.43, suggesting that VSGs absorb almost half of the UV energy that is absorbed by dust 
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Table 1. Total integration time and usage rate of short-exposure images in MIR 



Wavelength 


Total integration time of images 
sec 


used 


Percentage of images used 

% 


7.7 




2.7 


23 


8.8 




5.2 


43 


9.7 




4.9 


41 


10.4 




24.9 


42 


11.7 




53.7 


53 


12.3 




68.2 


47 




120 150 200 



T,o„t(MIR) (K) 

Fig. 7. — Ratio of emissivities measured in the MIR and in the i^band compared with color 
temperature measured in the MIR. Error bars are the la uncertainty of each measurement. 
The solid line shows the relation £eont(-^band)/£cont( MIR) = (7;ont(MIR)/651 K)^-^ (eq. [5]). 
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grains in thermal equilibrium. The uniformity of the ratio F(i^band)/F(MIR) suggests that 
dust in thermal equilibrium and VSGs are mixed well around the central engine. 

4.2. Spectral profile change of the silicate absorption feature 

At 10.4 /im, the silicate absorption towards the central engine is deeper than expected 
from the standard silicate feature. Figure 1 compares the observed SEDs with expected 
spectra from the IR excess towards /i Cep (Roche & Aitken 1984). There are a number of 
similar wavelength shift of the silicate feature in the literature. The interferomctric spectrum 
by JafFe et al. (2004) also shows a shift of the peak of the absorption feature towards longer 
wavelength. Silicate emission feature detected towards a type 1 AGN (Sturm et al. 2005) 
and quasars (Siebenmorgen et al. 2005) is shifted to a longer wavelength. 

The spatially integrated SED in Figure 1 is consistent with the spectrum observed 
through the 14"x20" aperture of the ISO-SWS (Sturm et al. 2000) which shows the normal 
spectral profile of the silicate absorption feature. This suggests that the spectral profile is 
normal as a whole and there is a spatial change of the profile. 

To measure the optical depth of the feature without being affected by the change, we 
fitted the MIR SED with a gray body emission model without using the 10.4 /im image 
(§2.3). If the 10.4 fim image were included in the grey body fitting, the resulting r9.7 would 
have been higher than expected from the literature (Roche et al. 1984; Tomono et al. 2001), 
especially near the central engine where T9.7 would have to be significantly adjusted for the 
low flux at 10.4 iim. 

Figure 8a shows spatial distribution of optical depth in the 10.4 fim filter rio.4. It is 
measured with comparing the 10.4 /xm image with the model continuum surface brightness 
fitted from the other images. The shape of the 10.4 iim absorption area is different from the 
9.7 iim absorption area (Figure 26). To illustrate the difference, we compared the measured 
T10.4 with that expected from the measured T9.7 and the IR excess towards /i Cep (Roche & 
Aitken 1984). Figure 8b shows the excess optical depth at 10.4 //m (AT10.4) near the central 
engine and to the north north east. 

The change of the spectral profile might be due to change of dust properties as suggested 
by Siebenmorgen et al. (2005) and Sturm et al. (2005) for type 1 AGNs and quasars. The 
area with large Ario.4 in Figure 86 coincides with the pie shaped high emissivity area (§ 4.4) 
seen in Figure 2d. This suggests that the spectral profile change of the sihcate absorption 
feature occurs in dense regions. 
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Aa (arcsec) 

Fig. 8. — Spatial distribution of the silicate absorption feature at 10.4 fim and 9.7 ^m. 
{a) Optical depth measured at 10.4 /im. Solid contours show rio.4 between 0.1 and 0.9 
(absorption; stronger with darker grey scale) with intervals of 0.1 while dashed contours 
show rio.4 between and —0.3 (emission; light). (6) Excess optical depth at 10.4 /xm (Ario.4; 
see text) for the jj, Cep extinction curve. Sohd contours show Atio.4 between 0.05 and 0.35 
(larger T10.4; darker grey scale) with intervals of 0.1 while dashed contours show Atio.4 below 
—0.05 with intervals of 0.1. Locations of the 5 GHz sources (GaUimore et al. 1996, open 
circles) and the [O 111] clouds (Evans et al. 1991, filled circles) are indicated. 
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4.3. Spatial variation of 9.7 /xm and 3.4 //m absorption features 

It is known that optical depth of the absorption feature of sihcate dust and that of 
carbonaceous dust are proportional towards Galactic sources. Roche & Aitken (1984) showed 
that Ay /tc^ j = 18.5 towards Wolf-Rayet stars and B supergiants in our Galaxy. Pendleton 
et al. (1994) measured optical depths of the carbonaceous dust and concluded that Av/tsa 
is ~ 150 for sources in the Galactic center or ~ 250 for sources in the local interstellar matter 
(ISM). Prom these results, ratio of the optical depths of the absorption features T9.7/T3.4 is 
~ 8.1 for the sources in the Galactic center or ^ 13.5 for the local ISM. Towards AGNs, 
T9.7/T3A is expected to be smaller because of temperature gradient in the dusty torus in 
the line of sight (Imanishi 2000). In fact, they measured T3 4 towards NGC 1068 in the 3'.'8 
aperture, compared it with Tgj measured in the 4'.'7 aperture by Roche et al. (1984), and 
yielded T'^.j/t^a = 4.3. With the data with higher spatial resolution, spatial variation of the 
ratio might show a trace of the dusty torus. 

Our results show significant spatial variation of the ratio rsiiicate/T3.4 (Pigure 9), with a 
complexity not expected to be produced by a dusty torus. As shown in Pigure 1, the 9.7 fim 
silicate feature is shifted towards longer wavelengths at some locations: optical depth of 
the feature is deeper at 10.4 //m than at 9.7 /xm (§ 4.2). To minimize influence from the 
wavelength shift, we defined TgiUcate ^ either T9.7 or T10.4 whichever has the higher absolute 
value. At the central peak, the ratio from our observation is between that measured by 
Imanishi (2000) and the Galactic vahics. In the slit C, the ratio increases from the central 
peak towards the north and approaches the Galactic values. In the eastern slits, the ratio 
is fiatter and closer to that measured by Imanishi (2000). In the south, the ratio is smaller. 
Pigure 4 shows that the silicate feature is observed as emission in the south. 

The emission of the silicate feature in the south of the central peak implies a temperature 
gradient in the line of sight: e.g., an optically thin layer of warm dust is present in front of 
cooler continuum source. On the other hand, in the north, the change of the ratio rsiiicatc/''3.4 
may be explained either with a temperature gradient in the line of sight: warmer dust clouds 
are located behind cooler dust clouds, or with differences in dust composition. Spatially 
resolved MIR spectra are needed to disentangle the possibihties. 

4.4. Pie shaped weirni dust cloud 

With the spatial resohition of 26 pc on the sky, the plausible dusty torus is not detected. 
Instead, we found a pie shaped area of high MIR emissivity in the north up to about 50 pc 
from the central engine (Pigure 2d). The pie shaped area is in the west of the radio jet. 
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Fig. 9. — Ratio of optical depths of the dust features. The numerator TsiUcate is either 
Tg 7 or rio.4, whichever had the higher absolute value. The vertical dotted lines from left 
to right show the ratios measured on the 3'.'8x3'.'8 aperture (Imanishi 2000), on Galactic 
center sources, and on local ISM (Pendleton et al. 1994), respectively. Error bars are the la- 
uncertainty of each measurement. The vertical line at top left shows the spatial resolution 
(FWHM) of the data. 
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There is no other hint of a channel supplying material to the central engine in the IR images 
in the literature and in this study. Therefore, we speculate that the area might be a channel 
that supplies material to the central engine. 

The lifetime of this cloud can be estimated assuming that material falls down into 
the central engine through the cloud. The dust mass of the cloud is estimated at around 

4.6 X 10^ Mq integrating e^ont{MlR) over the cloud where £cOTii(MIR) > 10 and dividing 
with Kabs (§ 2.3). Assuming the dust-to-gas ratio of 100 (Contini & Contini 2003), total 
mass of the cloud is ~ 4.6 x 10'^ Mq. Mass accretion rate to the central engine is estimated at 
0.05 Mq year~^ from total luminosity L = 2.2 x 10^^ Lq (Telesco & Harper 1980, at luminosity 
distance of 15 Mpc) and a typical value of mass-to-luminosity conversion efficiency of a black 
hole L ~ 0.3c^dM/dt (Eardley & Press 1975). The estimated total mass and accretion rate 
yields a hfetime of the cloud to be on the order of 10^ years. 

5. CONCLUSION 

Wc observed an area of 2'.'8 (200 pc) around the central peak of NGC 1068 in the MIR 
(8.8 - 12.3 /xm) and in the i^band (3.0 - 3.9 /xm). The shifted-and-added MIR images 
have a spatial resolution of 0'/37 or a projected distance of 26 pc while the L-band spectra 
are taken with a seeing hmited spatial resolution of O'.'S or 22 pc. Prom these data, we 
derived grey body parameters: color temperatures of the continua and emissivities, which 
are proportional to column densities of dust emitting the continua, at each spatial element. 
Moreover, optical depths of spectral features of silicate around 9.7 //m and carbonaceous 
dust around 3.4 /im were also derived (§§ 2.3, 3.3). 

The extended continua over 100 pc to the north and to the south of the central source 
are detected both in the MIR and in the L-band. We found that the MIR continuum is 
mainly emitted from dust in thermal equilibrium with radiation from the central engine 
while the L-band continuum is emitted by VSGs (§ 4.1). The observed SED in the MIR 
suggests that the peak wavelength of the 9.7 /xm silicate absorption feature is shifted to a 
longer wavelength at some locations (§ 4.2). The ratio of the optical depths of the silicate 
and carbonaceous dust features show complicated spatial distribution (§ 4.3). There is a 
pie shaped area of enhanced MIR emissivity extending about 50 pc to the north from the 
central engine. The morphology of the cloud leads us to speculate that the area is a channel 
that feeds material into the central engine (§ 4.4). 

We would like to acknowledge Miwa Goto for help in reducing the IRCS data. We 
would also like to thank the anonymous referee for the helpful comments which significantly 
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